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I examine some aspects of tlie interaction between the massive star Eta Cari- 
nae and its companion, in particular during the echpse-hke event, known as the 
^ spectroscopic event or the shell event. The spectroscopic event is thought to occur 

>- ■ when near periastron passages the stellar companion induces much higher mass 



loss rate from the primary star, and/or enters into a much denser environment 



cn 
(N 

^ ' around the primary star. I find that enhanced mass loss rate during periastron 

(^ ■ passages, if it occurs, might explain the high eccentricity of the system. However, 

^ ■ there is not yet a good model to explain the presumed enhanced mass loss rate 

■^ ' during periastron passages. In the region where the winds from the two stars 

Q^' collide, a dense slow flow is formed, such that large dust grains may be formed. 

Q ■ Unlike the case during the 19th century Great Eruption, the companion does 

-<— > ' not accrete mass during most of its orbital motion. However, near periastron 

ci ■ passages short accretion episodes may occur, which may lead to pulsed ejection 

of two jets by the companion. The companion may ionize a non- negligible re- 
^ ■ gion in its surrounding, resembling the situation in symbiotic systems. I discuss 

(Tj ■ the relation of some of these processes to other astrophysical objects, by that 

incorporating Eta Car to a large class of astrophysical bipolar nebulae. 
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1. INTRODUCTION 

1.1. The Binarity of r^ Carinae 

The similarity of the bipolar morphology of the Eta Carinae [r] Car) nebula— the Homunculus- 
(e.g., Ishibashi et al. 2003; Smith et al. 2004b) to bipolar nebulae observed in symbiotic neb- 
ulae (e.g., He 2-104 Corradi & Schwarz 1995; Corradi et al. 2001) and some planetary nebula 
(PN; e.g., NGC 3587 [PN G148.4+57.0] e.g., Guerrero et al. 2003) hints at a common process 
for the formation of the bipolar structure. Some similarities between the structure and spec- 
trum of 1] Car and PNs and symbiotic systems was noted before. For example. Smith (2003) 
points out to some similarities of the nuclear spectrum of the young bipolar PN Menzel 3 
[PN G331.7— 01.0)] with symbiotic systems and with rj Car, and to some common properties 
between the bipolar nebula of Menzel 3 and that of rj Car. Symbiotic systems are interact- 
ing binary systems, as was the progenitor of the bipolar PN NGC 2346 (PN G 215.6+03.6) 
which has a central binary system with an orbital period of 16 days (Bond 2000). This 
suggests that the binary companion to the primary of rj Car is responsible for the bipolar 
structure of its nebula (Soker 2001b; Soker 2004a,b). In these papers cited, I argue that 
the mechanism that most likely shaped the Homunculus, which was formed by the 20 years 
Great Eruption a century and a half ago (Davidson & Humphreys 1997), was two jets (or 
a collimated fast wind; CFW) blown by the companion during the Great Eruption (Soker 
2001b). The companion accretes the ejected mass from the mass-losing primary star. I also 
attribute the present fast polar wind found by Smith et al. (2003a) to interaction with the 
binary companion (Soker 2003). 

In those papers I also discuss problems with some single-star models for the shaping of 
the wind and circumstellar matter of r] Car (e.g., Langer, Garcia-Segura, & Mac Low 1999; 
Maeder & Desjacques 2001; Dwarkadas & Owocki 2002; Smith et al. 2003a; Gonzalez et al. 
2004; van Boekel et al. 2003). For example, Stothers (1999) found that the rotation does not 
much affect the instability of luminous blue variables, a view strengthened by Soker (2004a). 
In Soker (2004a) I discuss how the new finding of Smith et al. (2003b) of a more massive 
Homunculus further makes a single-star model for the formation of the bipolar nebulae of r] 
Car unlikely. After the paper on the problems that single star models have in explaining the 
bipolar structure of Eta Car was published (Soker 2004a), Matt & Balick (2004) proposed 
a model where the required rotation is somewhat slower than in most previous single-star 
models, ~ 0.2 times of the break-up rotation speed. For two reasons I consider their model 
to be problematic. In their model the progenitor of Eta Car has a huge magnetic field, i.e., 
the magnetic field pressure on the photosphere is larger than the thermal pressure; for that 
I consider their value of the magnetic field unrealistic. The second problem in their model 
is the rapid slow down of the envelope. In Matt & Balick (2004) model the magnetic field 
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keeps the wind in corotation with the envelope to several stellar radii. Hence, the specific 
angular momentum of the wind is very large, and rj Car rotation velocity will slow down on a 
short time. It is true that in the impulsive tidal interaction theory (rather than the quasi- 
static theory) the companion can spin-up the primary to a rotation speed of > 0.1 times its 
break-up rotation speed (Ivanov & Papaloizou 2004). However, during the Great Eruption 
the mass loss from the primary star had substantially spun-down the primary rotation. 

The binary model, of course, cannot account for all properties of 77 Car (Davidson 1999, 
2000), and the evolution and structure of the mass-losing primary star provide an explanation 
for many properties of rj Car. Such properties, i.e., evolutionary time scales and mass loss 
rate, should include stellar rotation as well (e.g., Heger & Langer 2000; Maeder, & Meynet 
2003; Meynet & Maeder 2003; Hirschi, Meynet, & Maeder 2004). 

The binary nature of rj Car is inferred from the so called spectroscopic event— the fading 
of high excitation lines (e.g., Damineli et al. 2000). For more observational support for the 
presence of a binary and its properties see, e.g., Damineli (1996), Damineli, Conti, & Lopes 
(1997), Damineli et al. (2000), Ishibashi et al. (1999), Corcoran et al. (2001a,b, 2004b), 
Pittard & Corcoran (2002), Duncan & White (2003), Fernandez Lajus et al. (2003), and 
Smith et al. (2004a). Since the spectroscopic event could be a result of mass-shell ejection 
by the primary star (Zanella, Wolf, & Stahl 1984; Davidson et al. 1999; Smith et al. 2003a; 
Martin & Koppelman 2004) it is also called a shell event. The main motivation to assume a 
shell ejection is in a single star model; in the binary model the shell ejection is not necessary, 
although it may occur. The periodicity of the spectroscopic event is seen in many wave bands, 
from the IR (e.g., Whitelock et al. 2004), to the X-ray (Corcoran et al. 2001a; Corcoran 2004; 
Corcoran et al. 2004a,b). The presently observed orbital period is ~ 5.5 yr (2023 ± 3 days 
as given by Whitelock et al. 2004). For primary and companion masses of Mi = I2OM0 
and M2 = 30Mq, respectively, the semi major axis of the orbit is a = 16.6 AU (there is 
not yet agreement on all the binary parameters, e.g., Ishibashi et al. 1999; Damineli et al. 
2000; Corcoran et al. 2001a, 2004b; Hillier et al. 2001; Pittard & Corcoran 2002; Smith et 
al. 2004a). 

The present study reports the possible effects of the stellar companion, mainly during 
the spectroscopic events, with the goal of comparing these effects with similar ones which 
could have occurred during the Great Eruption, and to strengthen the link between rj Car 
and similar bipolar nebulae with binary central stars. There is not yet agreement on the 
role the companion plays. Smith et al. (2004a) begin their paper by raising the question 
of what role a companion plays. As in my previous papers, I attribute most, or even all, 
of the shaping of the nebula around rj Car to the companion and its interaction with the 
primary mass losing star. The different sections of the paper discuss different processes and 
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their relevance both to the spectroscopic events and to the wind shaping during the Great 
Eruption. I summarize in section (8). 



1.2. Accretion States in i] Carinae 

Before turning to the new parts, I describe the three different accretion states of rj Car 
that I refer to in this and the previous papers. The basic flow structure is as follows. Both 
stars, 1] Car (the primary) and its companion, blow winds. The winds collide, and at one 
point the momentum flux of the two winds exactly balance each other, forming a stagnation 
point. The stagnation point is located close, but not exactly on, the line joining the centers of 
the two stars (the stagnation point is not exactly on the line because of the orbital motion). 
The shocked material of the primary star cools very fast (Pittard & Corcoran 2002; see eq. 1 
in Soker 2003), namely, before the mass moves to a large distance from the stagnation point. 
The surrounding pressure then compresses the cooled post-shock gas to high densities. Most 
of the mass is blown by the primary star, hence the state of the secondary's wind is less of 
an interest here. 

When blowing the wind, in both stars radiation pressure on the escaping gas overcomes 
gravitational attraction. However, this might not be the case with the dense gas near the 
stagnation point. Gravitational force on the dense and slowly moving (much below escape 
velocity) gas there might become large enough to accrete part of the mass back. Whether 
accretion occurs at all, and if it does, which of the two stars accrete most of the gas depends 
on the state of rj Car. I distinguish between three states. 

1. The Great Eruption. The accretion during the Great Eruption was studied by Soker 
(2001b). During the Great Eruption the mass loss rate by the primary was very high, 
and the stagnation point was within the companion's (secondary star) Bondi-Hoyle 
accretion radius (see discussion following equation [2] in the next section). The impli- 
cation is that along the entire orbit the companion steadily accreted mass with high 
specific angular momentum. Accretion disk was formed, according to the collimated 
fast wind (CFW) model, and two jets (or a CFW) were launched. These jets shaped 
the two lobes which are now observed as the Homunculus. The companion is expected, 
according to the model, to accrete ~ 50% of the mass blown by the primary, and blow 
~ 10 — 20% of the the accreted gas in the jets. The detection of jets in massive young 
stellar objects (YSOs) (Davis et al. 2004) show that jets can be blown by massive 
stars. 

2. Apastron passages in present r] Car. In Soker (2003) I proposed that during 
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present apastron passages the primary itself can accrete a fraction of ~ 5% of the mass 
blown over an entire orbit. This should not be confused with the accretion fraction 
of ~ 50% by the companion during the 20 years Great Eruption. The assumption is 
that part of the cool post-shock primary's wind material stays bound to the system. 
Near apastron passages the primary orbital motion is slow, and since the the system 
spends most of the time near apastron, it was claimed that the primary itself, due to 
its larger mass, will accreted back most of the bound mass. The high specific angular 
momentum of the accreted gas implies the formation of an accretion disk, which might 
blow a CFW. 

3. Periastron passages in present rj Car. This state was not discussed before, and 
it is a new idea presented in more detail in the next section. Basically, the stagnation 
point of the present wind of rj Car is outside the accretion radius of the companion, 
and the companion is not expected to accrete mass. However, if one or two of the 
following occurs, then accretion is plausible. First, if a factor of ~ 20 enhanced mass 
loss rate near periastron passages occurs (Corcoran et al. 2001a), then the stagnation 
point is about the same distance from the companion as the accretion radius (see 
next section). Second, when the companion is close to 77 Car, the wind speed may 
not reach its terminal speed yet, in particular if tidal effects of the companion on 
the acceleration zone of the wind from 77 Car are important. The denser 77 Car wind 
increases the likelihood of accretion. In any case, the accreted gas is expected to be in 
the form of dense blobs (cooled in the post-shock region), leading to sporadic accretion. 
This type of flow deserve a 3D gas-dynamical numerical simulations. 



2. MASS ACCRETION AND WINDS COLLISIONS 

The main difference between the Great Eruption and the present flow structure of 1] 
Car is the strength of the primary's wind (see previous subsection for the different accretion 
states). The argument goes as follows (Soker 2001b): The Bondi-Hoyle accretion radius of 
the companion, for accretion from the primary's wind is 

2GM2 f M2 \ f Vr \-2 , 

^a^— r^ = 0.2 — ^ — -^-— AU, 1 



r 



vl ■ V30M^/ VSOOkms 







where Vr is the relative velocity between the wind and the companion. As discussed in Soker 
(2001b), Vr = 500 km s~^, is a reasonable value including both the wind and orbital velocities. 
The distance D2 of the stagnation point of the colliding winds from the companion along 
the line between the stars, is given by equating the ram pressures of the two winds pv"^. For 
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spherically symmetric winds 

D, = r(3{l + (3)-'c:^r(3 = a( / - ) I 1^ [ , (2) 

Vl + ecos6'y yMiViJ 

where r is the orbital separation, 6 is the angular distance along the orbit (^ = at perias- 
tron), and P = [{M2V2) / {MiVi)Y^'^ . In the second equality I assumed P <^ 1. The stagnation 
point should be compared with the accretion radius. For the present winds' parameters 
of rj Car, with no periastron enhanced mass loss rate, taken from Corcoran et al. (2001a), 
vi = 500 km s-\ Ml = 10"^ Mq yr~\ V2 = 2000 km s-\ and M2 = IO^^Mq yT'\ I find 
f3 ~ 0.63, and D2 ~ 0.6 AU. The accretion radius is smaller than D2, hence no continuing 
accretion takes place. During the Great Eruption, when mass loss rate was much higher, 
such that at periastron j3 ~ 0.02, and D2 < Ra, accretion took place near periastron, and 
indeed along the entire orbit if the primary wind velocity was lower (Soker 2001b). 

If we consider a factor of 20 enhanced mass loss rate lasting 80 days after periastron 
(Corcoran et al. 2001a; but see section 6), then f3 = 0.14 and D2 — 0.2 AU. The accretion 
radius is comparable to D2. Since the shocked material at and near the stagnation point 
slows down and cools fast (Soker 2003), it is possible that part of it will be accreted by the 
companion during the ~ 80 days of enhanced mass loss rate. If the enhanced mass loss rate is 
due to lifting of gas from the acceleration zone (section 6), then the primary equatorial wind 
will be very slow. A slow equatorial flow, v < 100 km s~^ is observed indeed in the vicinity 
of f] Car (Zethson et al. 1999). This will not change much the Bondi-Hoyle accretion radius 
near periastron because of the high orbital speed, and equation (1) applies. The distance 
D2 will increase by a factor of ~ 2 due to the lower speed, and the accretion radius will be 
somewhat smaller than the stagnation distance D2. Overall, I expect that accretion will not 
occur during most of the periastron passage, but sporadic accretion of very cold parcels of 
gas (and dust) from the wind collision region (see section 3) might occur. 

If we consider a slow wind of f ~ 100 km s~^, then to obscure X-ray emission during the 
spectroscopic event, the mass loss rate not need increase by a factor of ~ 20 as required in 
the X-ray model of Corcoran et al. (2001). A modest increase of mass loss rate, by a factor of 
~ 3 — 5, is sufficient, and only in and near the equatorial plane where wind collision occurs. 
In that case, the wind interaction is more complicated, as the orbital speed is higher than 
the primary's wind speed, such that it changes the geometry. The stagnation point is still 
D2 ~ 0.2 AU, but its location is in front of the companion in its orbital motion, rather than 
facing the primary. The relative velocity between the companion and the wind is somewhat 
smaller, such that the Bondi-Hoyle accretion radius is somewhat larger. Sporadic accretion 
and disk formation may occur in this case as well. 
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As the dense post-shock material flows toward the secondary star three relevant forces 
are acting on it: gravity, wind's ram pressure, and radiation force. For the parameters of the 
companion star, wind's ram pressure is not much smaller, and probably somewhat larger, 
than radiation pressure (depending on the unknown secondary luminosity). It is therefore 
adequate to consider the force due to wind's ram pressure and gravity. I consider a post- 
shock spherical blob of mass nib, density pb, and radius r?,, located at a distance r from 
the secondary of mass M2. I also take the blob's density to be B times the wind's density 
Pb = Bpw2- The magnitude of the gravitational force is given by 

fa = ^2 = ^:^ Y^fe^P-2, (3) 

where in the second equality I substituted for the blob mass. The force due to the ram 
pressure of the wind is given by 

fw2 = Pw2vlTirl. (4) 

The ratio of these forces is 

Jw2 o r^V2 r r 

where in the second equality I took the wind speed to be of the order of the Keplerian speed 
on the surface of the secondary star, of radius R2. For the present system, the first factor 
on the far right hand side of the last equation is R2/r > R2/D2 ~ 0.03, where the minimum 
value is obtained for the distance of the stagnation point D2 ~ 0.5 AU. For pressure 
equilibrium between the ram pressure of the secondary wind having V2 = 2000 km s^^, and 
a cool blob at a temperature of Tb ~ 10^ — W^K (see next section) the compression factor 
is -B > 3 X 10^. The requirement for the blobs to be accreted, i.e., F > 1, therefore, is that 
its size be r^ > O.OOlr. This condition is likely to be met, as the region near the stagnation 
point from which mass is accreted is much larger; many instability modes can be developed 
with sizes larger than O.OOlr = O.OOI-D2 near the stagnation point. 

As stated, the sporadic accretion during the spectroscopic event may lead to disk forma- 
tion, and jets blown by the companion. The dense accreted blobs implies a short period of 
high accretion rate. Therefore, if jets are blown they will be of short duration, i.e., high mass 
loss rate into the jets. If the jets are narrow, then their momentum flux (momentum per 
unit time per unit area) be larger than the momentum flux of both the primary's wind and 
the secondary's wind. Hence, the jets will expand to large distances and to the Homunculus. 
I would like to speculate that these jets could explain the high velocity component observed 
in the X-ray emission of the iron K-shell line at 6.4 kev (Corcoran et al. 2004a). It is 
possible that the origin of the strong He II 4686 line observed by Steiner & Demineli (2004) 
to reach a velocity of —400 km s~^ near periastron passages, is also related to the accretion 
flow or the outflow. 



3. DUST FORMATION IN WINDS COLLISION 

After colliding with the secondary's wind, the wind gas of the primary star shocks to 
a temperature of ~ 3 x 10^ K (Corcoran et al. 2004b, and references therein). The cooling 
time from this temperature to a temperature of ~ lO'* K is much shorter than the flow time 
(Pittard & Corcoran 2002; Soker 2003). The flow time is deflned as 

rf = - = 3.5 f 4^1 { :' , ) " days, (6) 

^ vi \l AUy VsOOkms-i/ "^ ' ^ ' 

where Vi is the primary's wind speed, and Di is the distance of the stagnation point from 
the primary star. However, because the cooling time is much shorter than the flow time, the 
post-shock gas is not accelerated much, and the flow time out of the colliding region is much 
longer. For an outflow speed equals the speed of sound at ~ 10"^ K, the effective flow time is 
an order of magnitude longer rje ~ 1 month. For the wind parameters used in the previous 
section, e.g., an enhanced mass loss rate by a factor of ~ 20 during periastron passages 
(Corcoran et al. 2004), the hydrogen number density of the pre-shock wind at distance Di 
from the primary is 

-H^ = 2 X 10- { /^ \ { ;^ \ ' (^] " cm-. (7) 

\^10-4Mq yr-iy V500kms-i/ \l AU J ^' 

Behind the shock wave the gas cools and it is compressed by the ram pressure of the wind to 
much higher densities. Equating the thermal pressure of gas at 10^ K to the ram pressure of 
the wind, gives the density in the cool post-shock gas np(10^ K) ~ 2x IO^uhw ~ 4x 10^^ cm~^. 
Based on flgure 11 of Woitke, Kriiger & Sedlmayr (1996), I flnd that the cooling time of 
gas in the density range hh — 10^^ — 10^^ cm~^ from a temperature of 10^ K to 10^ K is 
< 1 month. This implies that as the post-shock gas flows away from the stagnation point 
(in the stage where no accretion to the companion occurs), it has time to cool to dust- 
forming temperatures. By then, the post-shock gas will be further compressed to densities 
of ~ 10^^ — 10^^ cm~^. These densities are much higher than in the gas where dust forms 
around AGB stars, resulting in more efficient formation of dust, e.g., large grain. Therefore, 
dust with properties different than the dust around AGB stars might be found in rj Car 
environment, in particular after the system emerge from a spectroscopic event. Based on 
infrared observations (e.g., Whitelock et al. 1994; Smith et al. 1998), and on extinction 
(e.g., Hillier et al. 2001) the presence of large grain in r] Car was proposed before (See Smith 
et al. 2003b for detailed discussion). 

Large dust grains are found in the disk around the binary system HD 44179, with one 
component being a post-AGB star, which is located at the center of the Red Rectangle, 
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a bipolar protoplanetary nebula. This binary system has an orbital period of Torb = 322 
days, a semimajor axis of asini = 0.32 AU, and an eccentricity of e = 0.34 (Waelkens et 
al. 1996; Waters et al. 1998; Men'shchikov et al. 2002). Large grains are also inferred 
in observations of other protoplanetary nebulae (e.g., Sahai, Sanchez Contreras, & Morris 
2004). It is commonly assumed that large dust grains are formed in a long-lived circumbinary 
disk around HD 22179 in the Red Rectangle (e.g.. Jura & Kahane 1999). In an earlier paper 
(Soker 2000b) I suggested that large grains can be formed in a slow and dense outflow, hence 
a long lived disk is not necessary (this view was strengthened by Men'shchikov et al. 2002). 
Men'shchikov et al. (2002) take for the grain-growth time to radius ag 

o-g \ ( Pg \ f nu \ -1 / T ^ '^'"^ 



^gg-0-2(^Q2cmv' V2gcm-3j Vl0i5cjn-3J VlOO kJ ^''' 

This shows that millimeter size grains can indeed formed near periastron passages of t] Car. 

4. INCREASING ECCENTRICITY 

The eccentricity e is reduced by tidal forces on a time scale called the circularization 
time and is defined as Tcirc = — e/e. Therefore, either tidal circularization is weak, or the 
high eccentricity of t] Car, e = 0.8 — 0.9, should be accounted for by a counter effect. In 
this section I study these effects in 7] Car. Here I take results from the standard theory 
of quasi static tides. However, the tidal interaction in high-eccentricity binary systems 
contains contribution from impulsive tidal as well (Ivanov & Papaloizou 2004). In a future 
work, devoted to tidal interaction, the tidal evolution of the r] Car binary system should be 
calculated with the impulsive tide contribution included. However, this is beyond the scope 
of the present paper. 

In the common tidal model in use, the equilibrium tide mechanism (Zahn 1977; 1989), 
the circularization time for evolved stars is derived by Verbunt & Phinney (1995). Here I 
use the equation as given by Soker (2000a), but change the scaling to fit rj Car. 
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env 



where L, R and Mi are the luminosity, radius, and total mass of the primary star, M^ 
is the primary's envelope mass, and /c ~ 1 is a dimensionless function of the eccentricity. 
The envelope is assumed to be convective. If it is radiative, then the circularization time is 
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longer. For a slow rotating star the function /c(e) is given by (Hut 1982) 

m ^ (1 - ef/^ (l + ^e^ + ^e^ + le^) , (10) 

where I moved a factor of (1 — e)^ from the usual definition of /(e) to the last term producing 
a more transparent equation, and when e = 0.9, fc = 0.17. This shows that tidal interaction 
will not reduce much the eccentricity, even for a more massive envelope, and even if during 
part of the time the primary's radius is larger than its present value. 

On the other hand, the eccentricity may increase because of the interaction, via enhanced 
mass loss rate at periastron passages or via a circumbinary disk. The change in eccentricity 
due to an isotropic mass loss (the derivation is applicable for an axisymmetric mass loss as 
well) is given by (Eggleton 2005) 

<5e = ^(e + cos^), (11) 

where SM is the mass lost from the binary in the stellar wind at the orbital phase 6 (hence 
SM < 0), M is the total mass of the binary system, and 6 is the orbital angle, measured from 
periastron, of the position vector from the center of mass to the secondary. The derivation 
of equation (11) assumes that SM{6) = SM{—6). As in Soker (2000a) I assume that in 
addition to its constant mass loss rate over the orbital motion M^, the primary star loses 
an extra mass 6Mp in a short time during the periastron passage, cos^ = 1. The total mass 
being lost in one orbital period Torb, is AMo = M^jTorh + SMp. I define the fraction of the 
mass being lost at periastron 

P-^^ (12) 

Under the assumption that the fraction of mass lost at periastron passage (3 does not change 
during the evolution, the relation between initial and final eccentricity is (Soker 2000a) 

where Cj and Mj are the initial eccentricity and total binary mass, respectively, and M is the 
final total mass. For example, if on average during the entire evolution (mainly before the 
Great Eruption), half of the mass was lost during periastron passages, i.e., f3 = 0.5, and the 
total system mass changed from Mi = ISOMq to a final (present) mass of M = 14OM0, then 
(1 + e)/(l + Cj) = 1.13. For example, when Cj = 0.7, e = 0.9. This shows that an enhanced 
periastron mass loss rate, if occurring, can increase somewhat the eccentricity. In any case. 
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equations (9) and (13) explain why the system eccentricity was not reduced substantially by 
tidal forces. 

Mass transfer can also change eccentricity. The change in eccentricity due to a mass 
(5Mtran transferred from the primary to the secondary is given by (Eggleton 2005) 

5e = 2<5Mtra„ f^ - ]i^) (^ + ^°'^)- (14) 

Since enhanced mass transfer may occur during periastron passage, we see that the eccen- 
tricity will decrease if Mi > M2, as is the case in r/ Car. Since for t] Car Mi ^ M2, we 
find from equations (11) and (14) that in order for the eccentricity not to decrease, the sec- 
ondary must not accrete more than a fraction of ~ M2/2M1 of the mass lost near periastron 
passages. During the Great Eruption, this situation may have been different because of the 
dense and somewhat slower wind (Soker 2001b), and the companion accreted mass along 
its entire orbit. This tends to reduce the orbital separation, but does not change much the 
eccentricity. 

As I mentioned in earlier papers (Soker 2003), t] Car is not unique. The most relevant 
system regarding this section is the binary system HD 44179 (period of Torb = 322 days, a 
semimajor axis of asini = 0.32 AU, and an eccentricity of e = 0.34; Men'shchikov et al. 
2002), which is located at the center of the Red Rectangle, a bipolar protoplanetary nebula 
(Cohen et al. 2004). In Soker (2000a) I argue that enhanced mass loss during periastron 
passages can account for the non-zero eccentricity of HD 44179. Waelkens et al. (1996), 
Waters et al. (1998), and Men'shchikov et al. (2002) attribute the non-zero eccentricity to 
the presence of a disk around the binary system. In this "external disk" mechanism, tidal 
interaction between the binary system and the circumbinary disk enhances the eccentricity, 
as is the model to explain eccentricities in young stellar binaries (Artymowicz et al. 1991; 
Artymowicz & Lubow 1994). 

Whether the mechanism for enhancing the eccentricity is an enhanced mass loss rate 
during periastron passages, or the external disk mechanism, the mechanism appears to be 
the same in r] Car and the Red Rectangle. The periastron mechanism implies that enhanced 
mass loss rate at periastron passage was a strong effect. The external disk mechanism implies 
that a relatively massive disk might be present, or might have been present before the 
Great Eruption, close to the binary system of t] Car. 

The time required for the secondary to spiral in to the primary is Tin = rcirc/16. This 
is a long time in 7] Car, during which the orbital separation increases because of mass loss. 
Hence, I do not expect the secondary to enter the primary envelope, unless the primary's 
envelope substantially swells. 
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TIDAL SPIN UP 



The synchronization time between orbital angular velocity and primary's spin angular 
velocity is r^yn ~ (1 + M2/Mi)(M2/Mi)-^(J/Mii?2)(i?/a)Vcir, where I is the primary's 
moment of inertia. For a crude estimate, I approximate the envelope density profile of 77 
Car by p oc r^^, as for stars on the upper AGB, where r is the radial distance from the 
star's center. For that density profile / = (2/9)Menv-R^- The scaled expression for the 
synchronization time reads 



'''° Vo-2/ Vio^^0/ Viooi?©; \\Mq) V0-25Mi 

For a slowly spinning star the function /s(e) reads (Hut 1982) 



a(l — e) 
3.6i? 



yr.(15) 



/.(e) - (1 + e)-^ (1 + f e^ + f e^ + ^e^) ; (16) 

when e = 0.9, /s(0.9) = 0.23. The short synchronization time implies that even without 
mass loss, which further slows down the envelope (Soker 2004a), the primary star will reach 
synchronization with the orbital motion during secular evolution. Namely, the primary will 
spin with a period of ~ 5 yr, which is ^ 10% of its break-up angular speed. The discussion 
above is based on the standard quasi-static tidal theory. However, in the impulsive tidal 
interaction theory, the companion can spin- up the primary to a rotation speed of > 0.1 
times its break-up speed (Ivanov & Papaloizou 2004). In any case, the mass loss from the 
primary star in the Great Eruption will substantially slow down the primary rotation, and 
the mass loss geometry is not expected to deviate much from spherical one (Soker 2004a). 
It is the binary companion which is behind the slow, dense equatorial mass flow, and it is the 
companion that ejected the collimated fast wind (jets) which inflated the two lobes during 
the Great Eruption (Soker 2001b). 



6. ENHANCING MASS LOSS RATE 

A suggestion that enhanced mass loss rate near periastron passages accounts for the 
spectroscopic event (Davidson 1999; Ishibashi et al. 1999; Corcoran et al. 2001) appears 
fundamental to some models as well as influencing some key issues in this paper. It is 
discussed in this section. 

What is the mechanism by which the companion increases, if it does, the mass loss 
rate near periastron passages? Some works related to other astrophysical objects propose 
enhanced mass loss rate as a result of the presence of a companion. The effect is basically 
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that the companion spins up the primary or rises a tidal bulge. A commonly used expression 
for the enhanced mass loss rate in binary systems is 






--(t 



;i7) 



where Mi, Ri, and Li, are the mass-losing giant stellar mass, radius, and luminosity, re- 
spectively, Rl is the radius of a sphere that has the same volume as the Roche lobe of the 
primary mass losing star, and Ai, Bl and 7 are constants. Different values for Bl and 7 are 
quoted in the literature: Tout & Eggleton (1988) use B^ = 10^, 7 = 6, and if Rg/Ri > 0.5 
then they set 0.5 for this ratio. Han et al. (1995) prefer B^ ^ 500 and 7 = 6, while Han 
(1998) argue for Bl ~ 1000 and 7 = 6. Frankowski & Tylenda (2001) argue for a more 
complicated expression, which basically has a very low value of Bl and 7 = 3; however, their 
numerical calculations yield a much faster increase in the mass loss rate as the giant is close 
to filling its Roche lobe. 

The derivation of equation (17), which uses the Roche lobe radius, assumes that the 
mass losing star corotates with the orbital motion. This cannot be the case in a highly 
eccentric orbit such as in the case of r] Car, where the orbital angular speed is much larger 
than the primary spin's angular velocity near periastron passages. Instead I take the point 
where the gravity forces of the two stars are equal. For an eccentricity of e = 0.9 and 
semimajor axis of a = 16.6 AU, at periastron passage r = 360Rq, and the gravity of the two 
stars are equal at a distance from the primary's center of ri ~ 240Rq. The radius of rj Car 
is -Rl ~ IOORq, such that Ri/ri ~ 0.4. Replacing Ri/Rl by Ri/ri in equation (17), gives 
the enhanced mass loss rate values of ~ 40, 3, and 5, for the parameters of Tout & Eggleton 
(1988), Han et al. (1995), and Han (1998), respectively. The prescription of Frankowski & 
Tylenda (2001) gives negligible enhanced mass loss rate. 

Frankowski & Tylenda (2001) try to find a mechanism to enhance the mass loss rate 
and find that it is not clear how the mechanism can substantially enhance the mass loss rate 
for Rl > 2Ri. One possible way out of this problem is by using a mechanism that does not 
affect the stellar surface as much as it affects the region from which the wind is accelerated. 
This region in giants can be quite large, in particular if dust formation occurs. Shocks from 
stellar pulsation can enhance or trigger dust formation (Woitke, Goeres, & Sedlmayr 1996). 
It might be, therefore, that the role of the companion in the case of rj Car is to make dust 
formation favorable above the atmosphere of rj Car, by enhancing gas density or inducing 
shocks. In the case of enhancing density, the relevant radius for the tidal interaction is not 
-Rtid = -Rl) but rather the radius where dust might form, which is -Rtid ~ 1.5 — 2.5-Ri. Since 
the height of a tidal bulge is proportional to (i?tid/''")^, where r is the orbital separation, an 
increase by a factor of ~ 2 in i?tid can make a huge difference. Is it possible that dust can 
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form as close to the surface of a hot star as the primary star of rj Car is (~ 20, 000 K)? I note 
that dust formation occurs near the surface of hot R Coronae Boreahs (RGB) stars. RGB are 
rare hydrogen-deficient carbon-rich supergiants which undergo very spectacular declines in 
brightness of up to 8 mag at irregular intervals as dust forms along the line of sight (Glayton 
1996). The hot RGB stars have an effective temperature of ~ 18,000 K, a luminosity of 
L ~ 10'^Lq, hence a radius of i? ~ IORq. Their surface temperature and gravity are similar 
to that of 1] Gar, but they are carbon rich, facilitating a lot dust formation. In any case, hot 
RGB stars show that dust can be formed close to the surface of hot stars. 

There is another difficulty. Even if the wind is slow, Vi ~ 100 km s~^, the wind 
momentum flux to radiation momentum flux ratio for an enhanced mass loss rate of 2 x 
IQ-^Mq yr-^ and luminosity of Li = 3 x IO^Lq, is quite high, C = MiVi/{Li/c) ~ 3. That 
is, the momentum flux in the wind is quite high, and it is not easy to explain this high 
momentum by radiation pressure. As stated in Section 2, the enhanced mass loss rate for a 
slow wind might be larger by only a factor of ~ 4, and near the equatorial plane. If the mass 
loss is spherical, then C < 1- I^i these circumstances, the supposition is that the enhanced 
mass loss rate occurs near the equatorial plane, where less radiation is available for mass 
ejection. 

Whitelock et al. (2004) argue that dust absorption has negligible role in the spectroscopic 
events, hence dust cannot absorb and scatter a large fraction of the radiation emitted by the 
primary star. 

Gonsidering the problems mentioned above, it is quite possible that the spectroscopic 
event is mainly a coverage of the region around the companion, as proposed by Whitelock et 
al. (2004). As discussed in Section 2, large quantities of dust could be formed in the region 
near the companion, possibly engulfing the companion as it deflects the flow around itself, or 
even accreting mass. Most of the radiation from the primary star is not blocked, but some 
fraction of the radiation is blocked. To estimate the magnitude of these effects, in the next 
section the ionization of the gas by the companion is discussed. 

This section is summarized by emphasizing again that it is not clear whether substantial 
increase in the mass loss rate from the primary occurs near periastron passages, or whether 
the secondary enters a more or less static dense medium. If the mass loss does increase, 
there is a larger question regarding the mechanism behind the increased mass loss rate. This 
question is an important open question regarding other astrophysical objects, for example, 
the Red Rectangle. 
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7. IONIZATION BY THE COMPANION 

The ionization of gas in the nebula of rj Car was discussed in a previous paper (Soker 
2001a). In that paper I propose ionization shadows as an explanation for the formation of 
the long and narrow strings of rj Carinae (see Weis, Duschl & Chu 1999), and I show that the 
companion could play a significant role in the ionization of the strings' surroundings. The 
main ionization was attributed to the time the companion is away from periastron, which is 
most of the orbital period. Here I examine the situation near periastron. 

The flow structure is very complicated near periastron passages, and 3D numerical 
simulations are required to find the shape of the ionized region. I therefore use a simple 
flow structure (as in Soker 2001a), wherein the primary's wind flows undisturbed around the 
companion, and I consider only the ionization of the primary wind, neglecting the companion 
wind. Although ideal, the derived expression gives a good indication of the main effect. 
Because of the complicated flow structure near the equatorial plane, I examine the ionization 
in the polar direction from the companion. Let r be the distance between the companion 
and the primary, and h the distance above the poles of the companion, i.e., perpendicular 
to the equatorial plane. The density of the primary's wind above the companion is p{h) = 
Mi/[47r(r^ + h'^)vi]. The total recombination rate, assuming spherical geometry for the sake 
of the calculation, along the direction h is given by 

N= A-Kane{h)ni{h)h^dh. (18) 

where n, and Ug is the ion and electron density, respectively, and a the recombination 
coefficient. Performing the integration from inner distance h^ to the edge of the ionization 
front hi gives 

iV = 47ra^^(^l ^[H,{l + Hl)-^-H,{l + Hl)-' + i^n-^H,-i^n-'H,], (19) 
p p \ 47rf 1 / 2r '- -' 

where H = h/r. 

To find the distance of the ionization front in units of the orbital separation Hi, I as- 
sume that Hq <^ Hi <^ 1, which is correct near periastron when the companion enters a 
dense medium (see below). Equating A^ to the ionizing photon luminosity of the compan- 
ion 5*2, scaled as in Soker (2001a), and using the approximation above, yields, after some 
manipulation 

—2/3 
hi / r \ 1/3 / Vi \ 2/3 / Ml \ / ^2 



r ■ V2AU/ VlOOkms-i/ \ 2 x IO'^Mq yr^M \ 2 x lO^s s-^ 
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The numerical value in the last equation will not change if we take vi = 500 km s ^ (no 
slower wind near periastron), and at the same time increase the mass loss rate to Mi = 

Although the last equation was derived under strong simplifications and assumptions, 
it does provide an insight into the ionization zone. As long as the companion is away from 
periastron, and the wind parameters are the present ones, i.e., mass loss rate is Mi ~ 
1O~^M0 yr~^ and Vi ~ 500 km s~^, the ionization zone is large, hi > r, and an accurate 
solution of equation (19) is required. If, however, near periastron the companion enters into 
a much higher density zone, caused by enhanced mass loss rate or slower velocity, then the 
ionization zone shrinks to a size hi <^ r. This implies that the ionization zone could be 
hidden by dense primary's wind, in particular by the dense gas in the region where the two 
winds interact. 

Finally, I note that a fiow structure where a companion ionizes the primary's wind is not 
unique to 7] Car. For example, it occurs in symbiotic binary systems and some protoplanetary 
nebulae. In the protoplanetary nebula M2-9 there is a positional shift of bright knots in the 
inner nebular lobes on a period of 120 years (Doyle et al. 2000). This side to side departure 
from axisymmetry is explained in terms of a revolving ionizing source, i.e., an ionization 
by one star in a binary system (Livio & Soker 2001). Livio & Soker (2001) show that the 
interaction between the slow, AGB star's wind, and a collimated fast wind from the white 
dwarf companion clear a path for the ionizing radiation in one direction, while the radiation 
is attenuated in others. Similar fiow and radiation structure (but not exactly identical in the 
geometry and winds parameters) was qualitatively applied to the side to side variation in the 
UV images before and after the spectroscopic event in r] Car (Smith et al. 2004a). In rj Car 
the UV images are of the near UV, with energy below the hydrogen ionization threshold. 
This UV radiation can expand to larger distances than the ionizing radiation studied in this 
section. It is interesting to note that the general bipolar structure of M2-9 was compared 
to the bipolar structure of r] Car by Smith (2003), and the nuclear spectrum of M2-9 was 
compared to that of r^ Car by Balick (1989) and Allen & Swings (1972). 



8. SUMMARY 

The main goal of this paper is to further explore the role of the binary companion in rj 
Car, and by that to strengthen the link between rj Car and other binary systems losing mass 
at high rates and having bipolar nebulae around them. This, I hope, will motivate more 
3D gasdynamical numerical simulations of the interaction of the two stars and their winds. 
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A unified model for sliaping bipolar nebulae (nebulae composed of two lobes with an 
equatorial waist between them) around mass losing stars has been emerge in recent years 
(Soker 2004b). In this model (Morris 1987, 1990), the two opposite lobes are shaped by 
a binary companion accreting mass from the mass-losing primary, and then blowing two 
jets. The two jets then inflate the bubbles, as in many other astrophysical objects, such 
as clusters of galaxies (Soker 2004b and references therein). Note that I refer here only to 
nebulae having two large lobes. Most elliptical planetary nebulae which have no large lobes, 
for example, were not shaped by this process. A second goal was to compare the present 
effects of certain processes with the effects of these processes during the Great Eruption of 
a century and a half ago. 

The main finding of this paper can be summarized as follows. 

1. As far as the binary interaction goes, the main question regarding the spectroscopic 
event is on the effect of the companion on the mass loss process near periastron pas- 
sages. Does the companion increase the mass loss rate from the primary? If it does, 
by how much? How does it change the geometry of the mass loss process, i.e., does the 
wind becomes slower near the equatorial plane? As discussed in Section (6), periastron 
enhanced mass loss rate was discussed in relation to other astrophysical objects, e.g., 
in the central binary system of the bipolar protoplanetary nebulae the Red Rectangle. 
However, if mass loss rate is indeed enhanced, a mechanism is yet to be identified for 
that to occur. A promising mechanism could be one in which the companion tidaly 
facilitates dust formation (Section 6). 

2. One of the results of periastron enhanced mass loss rate is a departure from axisymme- 
try in the equatorial plane. In the first type of departure from axisymmetry the mass 
loss rate near periastron and apastron causes a permanent differences between the two 
sides of the nebula. This is observed in rj Car (e.g.. Smith 2002; see discussion in 
Soker 2001b), and in the inner region of the Red Rectangle (e.g., Tuthill et al. 2002; 
Miyata et al. 2004), which is known to harbor a binary stellar system at its center. If 
the major axis of the eccentric orbit is perpendicular to our line of sight (Smith et al. 
2004a), then the observed departure in rj Car is small. 

3. There is another type of departure from axisymmetry, where there is a side to side 
variation as the two stars orbit the center of mass. This is observed in rj Car in the UV 
(Smith et al. 2004a) in a manner resembling that of the protoplanetary nebulae M2-9, 
where the orbital period is 120 years (Doyle et al. 2000). The explanation for this side 
to side variation during the orbital motion is the opacity of the dense wind blown by 
the mass-losing star, which allows the radiation from the companion to propagate only 
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in the other direction. This radiation conies from a hot white dwarf in M2-9 (Livio & 
Soker 2001), and from an O star in rj Car (e.g., Smith et ah 2004a). 

4. The winds from the two stars colhde, and the post-shock primary's wind cools on a 
time scale much shorter than the flow time. Even without enhanced or slow wind, near 
periastron passages a dense, cold region is formed. Large dust grains, of the size of a 
millimeter size, could be formed there (Section 3). Large dust grains are observed in 
the Red Rectangle; the popular dust formation site there is a long-lived circumbinary 
disk, although a slow, dense flow as the site for large grain formation was suggested in 
the Red Rectangle as well. 

5. If mass loss rate is indeed enhanced near periastron passages, or the wind is slower, 
the companion could sporadically accrete mass from the dense, cold region in the 
interaction region of the two winds (Section 2). Such accretion events may lead to 
the formation of an accretion disk around the companion, which could then blow two 
jets. The jets might expand to large distances, revealing themselves as fast velocity 
components. In Section (2) I speculate that these jets may be connected to the high 
velocity component observed in the X-ray emission of the iron K-shell line at 6.4 kev 
(Corcoran et al. 2004a). I note that during the Great Eruption, the companion accreted 
mass along its entire orbit and probably blew the jets that shaped the Homunculus 
(Soker 2001b; see section 1.2 here). 

6. Periastron enhanced mass loss rate increases eccentricity. If, during the evolution, an 
average mass fraction of f3 was lost near periastron passages, then the eccentricity was 
increased according to equation (13). In any case, tidal interaction did not much de- 
crease the eccentricity (Section 4), and if the eccentricity was large initially, it remained 
so. High eccentricity, although not as high as in rj Car, is found also in the bipolar 
nebula the Red Rectangle. 

7. Because of short synchronization time (equation 15), the the primary spin period is 
< 5 yr, meaning that its rotation speed is much below its break-up rotation. Hence, 
effects of its rotation on the axisymmetrical mass loss geometry are expected to be 
small (Soker 2004a) 

8. During most of the orbit the companion is likely to ionize a large volume of the pri- 
mary's wind (Soker 2001a). This process of ionization by a companion occurs also in 
symbiotic binary systems. Some similarities of the nuclear emission of rj Car with 
symbiotic systems and planetary nebulae was noted before (e.g., Balick 1989; Smith 
2003). However, near periastron passages, if the r] Car primary's wind is slow and/or 
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mass loss rate is higher, the ionized region shrinks, such that it can be echpsed, par- 
tially or fully, by the primary's dense wind (Section 7), as the X-ray emission is eclipsed 
(Corcoran et al. 2004b). 
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